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We consider a simple model of a multidomain superconductor-ferromagnet-superconductor (SFS) 
Josephson junction. Sign-alternating magnetization M in domains leads to a spatial modulation of 
the phase difference 0(a;). Due to this modulation the Josephson critical current Ic may have a 
different sign depending on the ratio of the magnetic flux in a domain, 4:nAIa{2dF), to the magnetic 
flux quantum. This phase modulation, but not a nonmonotonic dependence of the local critical 
current density jc, may be the reason for oscillations of the current Ic as a function of the F layer 
thickness 2dF or temperature, observed in experiments. 



I. INTRODUCTION 



New states have been and observed in Josephson junctions in the last years. These states are characterized by a 
negative Josephson energy Ej — {Ich/2e){l — cos0), that is, under certain conditions the Josephson critical current 
Ic changes its sign becoming negative. This means that the ground state corresponds to the phase difference equal 
to TT (TT-state), but not to zero as it takes place in ordinary Josephson junctions. Such states have been observed in 
Josephson junctions of different types:l) in multiterminal SNS Josephson junctions [1], 2) in junctions consisting of 
two d-wave superconductors (see references in the review [2]), 3) in SFS junctions [3-5] (see also the review [7]), where 
S,N and F stand for a superconductor, a normal metal and a ferromagnet, respectively. The 7r-state in SNS junctions 
is created by passing a dissipative current through the N layer. This current leads to a nonequilibrium distribution 
function of quasiparticles in the N wire with respect to the equilibrium distribution function in the superconductors. 
In the d-wave superconductors (high Tc-superconductors) the sign of the order parameter A depends on the direction 
in space with respect to crystallographic axes. Therefore if two S/N (or S/I, where I is an isolator) interfaces have 
different, properly chosen orientations, then the critical current Ic, which is proportional to the product A1A2, may 
be negative. This occurs provided that the order parameter in one superconductor (Ai or A2) is negative. In 
SFS junctions the critical current Ic may be negative because the condensate (Gor'kov's) Green's function, which 
determines the current Ic (or to be more exact, the critical current density jc), oscillates in space changing sign. The 
sign reversal of jc in SFS junctions has been predicted a long time ago by Bulaevskii, Kuzii and Sobyanin [8], who 
considered electron tunneling between two superconductors via a magnetic impurity. Later this effect was studied in 
SFS junctions by Buzdin, Bulaevskii and Panyukov [9] (references to other theoretical papers on this subject is given 
in the review [7]). It was shown that the current density jc decays with increasing the thickness of the F layer 2dF and 
changes sign (damped oscillations of jc)- Such a behaviour of the critical current has been observed experimentally in 
Refs. [3-6]; the current Ic decays with increasing the thickness 2dF or temperature T in a nonmonotonic way changing 
sign. In the recent paper [6] a spontaniously circulating current in a superconducting ring with a SFS 7r-junction has 
been observed. 

In this paper we show that the damped Ic oscillations are not necessarily related to such a dependence of the local 
critical current density jc as it was anticipated previously. The total Josephson current Ij, which is measured in 
experiments, is an integral from the local current density jj — jcSm(j) over the whole area of a SFS junction (we 
choose a simple, sinusoidal form of the dependence jj{4>), but the conclusions we make are valid qualitatively in 
a general case). It is important to have in mind that the phase difference varies in space in the presence of a 
magnetic field, and in multidomain SFS junctions a spatial dependence of 4> arises even in the absence of an external 
magnetic field H^xt- It was already mentioned in Ref. [3,4] that the maximum value of Ic corresponds to zero external 
magnetic field Hcxt- One can assume that this may be related to a multidomain structure of the F film. Otherwise 
the maximum value of Ic would be shifted by a certain value of Hcxt for which the induction Hcxt + 47rA/ is zero, 
where M is the magnetization in a one-domain F film. We consider a simple model of a multidomain structure of the 
F film and show that even if the local current density jc is always positive, the critical current Ic changes sign when 
the in- plane magnetic flux in a domain <^a — 4:TrMQ{2dpa) is an integer of the flux quantum $0, where a is the domain 
width. The domain width a depends both on the thickness dp and temperature T if the screening of stray magnetic 
fields by the Meissner currents is taken into account [10-12]. Therefore the magnetic flux is changed with increasing 
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dp and temperature leading to damped oscillations of the critical current Ic- It will be shown that the dependence 

Ic{^a/^o) may be described by a Fraunhofer-likc pattern 



i"c(«'a/$o) = /cO sin(7r$„/$o)/(7r$a/$o) 



(1) 



where Ico = jcLx is the critical current for an uniform (one-domain) SFS junction per unit length in the y-direction 
(we assume that all quantities depend only on x and z; see Fig.l). This dependence describes the critical current for 
the case of domains with an alternating magnetization which takes constant values ±Mo in domains. 



II. MODEL AND BASIC EQUATIONS 



Consider a SFS Josephson junction in which the F layer consists of stripe domains parallel to the y-axis (see Fig.l). 
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FIG. 1. Schematic picture of the considered SFS junction with magnetic domains (la). Fig. lb shows the supercurronts (solid 
lines with arrows) in the upper superconductor screening the normal component of a magnetic field near domain walls. These 
currents in the lower superconductor (not shown) have the same direction. The domain boundaries are represented by the 
dashed lines. The dashed linos with arrows illustrate the stray magnetic field related to rotation of the magnetization M in 
the domain walls. The cross and dot in the circles show the directions of the magnetization vector in neighboring domains. In 
Fig.lc the solid lines with arrows show the screening supercurrents in presence of the Josephson coupling. 

To be specific, we consider the Bloch domains, that is, the magnetization vector M lies in the y — z plane changing its 
direction in domain walls of the width w. In the case of a narrow domain width the magnetization vector M is parallel 
to the y-axis, constant in each domain and rotates in the domain wall. The magnitude \M\ is assumed to be constant 
everywhere in the F film. The orientation of the magnetization M in the F film (in-plane or out-of-plane) depends 
on many factors such as the Curie temperature, the constant of magnetic anisotropy, the film thickness etc (see for 
example, [13,14] and references therein). In ultrathin F films (a few atomic layers) a transition between the in-plane to 
out-of-plane M orientation may occur [15]. We consider much more thicker F films, used in experiments [3-6], where 
the magnetization orientation is determined by the direction of easy magnetization. Therefore it is assumed that 
the axis of easy magnetization is parallel to the plane of the F film. Even in the absence of the Josephson coupling 
between the superconductors S, the Meissner currents along the y-direction are induced by the stray magnetic field 
H = {Hx, 0, Hz) (dashed lines with arrows in Fig. la). The magnetic field components H^^z can be easily found from 
equations 

d''HxAkn,z)/dz^ - KlH,,z{kn,z) =0 (2) 

where Hx^z{kn,z) = JI^^{dx/2a)Hx,z{x, z) exp{iknx) is the Fourier component of Hx^z{x,z), kn = 'Kn/a, n = 
0,±1,...;k^ = fc^ -h kI, = Xl is the London penetration depth (in the ferromagnet the penetration depth 
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^L.F may be taken infinite because the amphtudc of the condensate in the F layer is small and therefore the screening 
is weak). This equation is supplemented by the boundary conditions [16] 

=0, [H,]=^TiM,{K,z) (3) 

where the square brackets mean a difi'erence: [Hx] = Hxikn, dp + Q) — Hxikn, dp — 0). The component is not zero 
in the domain walls. This implies that the in-plane component of H is continuous across the S/F boundary and the 

normal component of H, which exists near the domain walls, experiences a jump. One can easily solve Eqs.(2) and 
find the fields Hx,z which are connected with the vector potential A : Hx = —dAy/dz, Hz = dAy/dx. The expression 
for Ay in the superconductor is 

Ay{kn, z) = -'ii^Mz{kn){fnikn)~^ exp(-K„(2; - dp)); z> dp (4) 

where /« = 1 + Kn/{knta,iih.0n),9n = kndp. The vector potential Ay in the lower superconductor (z < dp) has the 
same sign. This means that the stray fields Hx,z lead to the screening currents 

jy{kn,z) = -{c/ATT)K\Ay{kn,z) (5) 

which have the same direction in the upper and lower superconductors, but the opposite directions in neighboring 
domains (sec Fig. la). Because the Meissner currents flow in the same direction in both superconducting electrodes, 
they do not lead to a phase difference between the superconductors and do not inluence the Josephson current 
essentially. These currents may locally reduce the amplitude of the order parameter and therefore decrease the 
critical current Ic if the magnetization Af^ is strong enough. We will not discuss this simple effect. 

If the Josephson coupling between the superconductors is negligible (this case was considered in Ref. [10]), the 
magnetic field has only the components Hx^z which are determined by the vector potential Ay{kn, z). The component 
Hy is zero. However the components Ax,z{kn, z) of the vector potential are not zero. For example in one-domain case 

Ax{z) = AirMyZ in F (6) 
Axiz) = ±A-KMydp in S (7) 

In the absence of the Josephson current the components Ax.z in a multidomain SFS structure are found from the 
equations: H = V x A = and V • A = 0. However in this case the presence of components Ax,z does not lead to the 
superconducting currents and therefore to the appearance of an additional magnetic field H. The term proportional 
to Ax in the expression for the supercurrent is compensated by the term proportional to gradient of the phase (see 
Eq. (9)). If the Josephson current jj is not zero, additional screening currents jx,z and therefore the component Hy 
arise in the system (see Fig. lb). The component Ax affects the phase difference and therefore the critical current Ic 
in the junction under consideration. 

In order to calculate the critical current Ic, one needs to derive an equation governing the phase difference ^ in a 
multidomain SFS Josephson junction. For simplicity we assume that the thickness of the F layer is small: 2dp << a, w, 
where 2a is the period of the domain structure and w is the width of the domain wall (one can analyze a more general 
case, but the calculations become more combersome). We need an equation for the component Hy which is related 
to the local Josephson current density 

(V X H)z = dHy/dx\z=o = (47r/c)ie sin (8) 

We assumed the simplest form of the relationship between the Josephson current density jj and the phase difference 
(f), but this assumption is not essential. We also dropped a contribution to the gauge-invariant "phase difference" 
which stems from the vector potential and has the form /_^^ Azdz. One can easily show that this part is smaller 
than (f) by the parameter {a/dp) (we choose a gauge in which V • A = 0). Let's write down the a;-component of one 
of the Maxwell equations for the current density in the superconductor at z = ±dp 

-dHy/dz = {A7r/c)jx = kH-Ax + {^o/2n)dx/dx) (9) 

where k| is defined in Eqs.(2) and x is the phase of the order parameter. Subtract Eqs.(9) taken at z = +dp and 
z = —dp from one another (a similar method was used by one of the authors in Ref. [18] in the study of collective 
modes in layered superconductors), we get 

- [dHy/dz] = kI{~[Ax] + {^o/2Tr)d(P/dx) (10) 

where the square brackets means, as before, a jump across the F layer and <p = xi.d'p)—x{~dp) is the phase difference. 
The jump [Ax] is found from the equation 
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AiTMy = dA^/dz-dA^/dx (11) 

and is equal to: [A2;(A:„)] — 4:nMy{kn)dF in accordance with Eq.(6)(see Footnote [19]). The contribution of the second 
term is smaller by the parameter (di?/a). The component of the magnetic field Hy can be found from an equation 
similar to Eqs.(2). With account for the boundary condition (10) this equation acquires the form 

d'^Hyikn, z)/dz^ - KlHy{ka,z) = ,5(z)4 ( + ($o/27r)z/c„(/.fe) (12) 

Solving this equation for Hy and substituting the solution into Eq.(8), we obtain an equation for Fourier components 
of the phase difference (pk = <Pikn,z) 

+ {ku/kl) {sin (j){x))k = iA;„47rMj,(fc„)(2dF)(27r/$o) (13) 

where kj = -^/T&Tr^j^JTc^i^ is the inverse Josephson length. In the coordinate representation Eq.(13) has the form 

-d'^(p/dx^ + j dxiK{x - xi)sin0(xi) = -4,iTdF{2TT/^Q)dMy{x)/dx (14) 
where the kernel K{x — xi) is defined as follows 

K{x - x\) = (l/2a) ^(kk/kl) exp(-ifc„(a; - xi)) (15) 

n 

Eq.(14) describes the dc Josephson effect in a simple model of multidomain SFS junctions with a thin F layer. 



III. TWO TYPES OF DOMAIN STRUCTURES 



If the London penetration depth is small compared to the domain size a and the width of the domain wall ( 
Al << a,w ), then Eq.(14) is simplified. This condition means that the characteristic kn values are much smaller 
than kl- In this case K{x — x\) 5{x — xi) and Eq.(14) acquires the form 

-d'^cp/dx'^ + K^jsm<j){x) = -ATT{2d,F){2TT/<i>o)dMy{x)/dx (16) 

This equation differs from the standard Josephson equation only by the term on the right-hand side. One can study 
various properties of the SFS junctions described by Eq.(16) or by Eq.(14), but in this paper we analyze only the 
critical current and its dependens on different parameters (external magnetic field, the thichkness dp etc). First we 
consider the case of a periodic My{x) dependence. If the period 2a of this dependence is much less than the long 
Josephson length Kj^, then a solution for Eq.(14) (the relation between Al and a,w may be arbitrary) is 

</' = 0M + 0O, <pM = {2TT/<^o)M^dF) J dXiMy{xi) (17) 

where (t>M{x) is a function fast varying in space, (/>o is a constant (or a function smoothly varying over the period a). 
The total Josephson current (per unit length in y-direction) is 

Ij = jc / dxi s\u{(pQ + (Pm{xi)) (18) 
Jo 

where (pM is given by Eq.(17). Note that the weak Josephson coupling does not affect the domain structure, and 
therefore this structure can studied in the absence of the Josephson efi'ect. As we noted, the domain structure was 
analyzed theoretically in Ref. [10] for arbitrary dp and in Refs. [11,12] for thick F layers {dp » a)- It was shown 
in Ref. [10] that the period a depends on rf^ in a nonmonotonic way and the width of the domain walls w may be 
much less than or comparable with the domain width a. One has a step-like structure My (x) in the first case and an 
oscillatory structure in the second case. Consider two limiting cases. 
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Ic/Ico 




FIG. 2. Normalized critical Josephson current as a function of the normalized magnetic flux in a domain m = 7r$a(a)/$o for 
different parameters 7 = a/(5, where a is the averaged domain size and 5/\/2 is the dispersion of the domain size fluctuations. 



a) Step-like domain structure. The magnetization vector equals:M = (0, Mj,(a;), 0), where My{x) — -\-Mq for 
< a: < a, and My{x) — —Mq for ~a < x < 0. Outside this interval the dependence My{x) is periodically repeated. 
For this structure the Josephson current Ij is described by the formula 



Ico sin 410 



sin(7r^a/$o) 

7r$a/$0 



(19) 



where Ico = jcLx, *&a = 47rMo(2c?i?a) is the in-plane magnetic flux in one domain. Therefore the critical current Ic, 
which is given by Eq.(l) oscillates and decays with increasing This implies that the amplitude of the critical 
current oscillations decreases with increasing dp ov temperature T because the period of the domain structure 2a 
depends on Xl{T) (see the theoretical papers [10-12] and the experimental paper [17], where it was shown that the 
domain structure is changed with changing T) . 

b) Oscillatory domain structure: My[x) — Afosin(fcox), fco = T^/a. 

In this case the critical current is equal to 



IcQjQ{4>M{a)) 



(20) 



where Jo is the Bessel function of the zeroth order. In both cases the behaviour of the critical current Ic as a function 
of $a is qualitatively the same: the current Ic decreases with increasing cpMio-) and changes sign. 

In our model of a periodic domain structure the action of an in-plane external magnetic field Hext on Ic can be 
easily analyzed. In the presence of the field Hext the phase difference equals (j){x) — (f>o + (j)H{x) + (t)M (x) , where (j)M{x) 
is given by Eq.(17), 0o is a constant and (/)_f/(x) = 2XLxHcxt{'^T^/^o)- If the domain size a is much less than the length 
of the junction Lx, the averaging over the period of the structure can be done as before at a fixed coordinate x, and 
we arrive at Eq.(19) in which one has to replace 4>o fpa + 4'h{x)- The final averaging over Lx yields, for example, 
in the model of a step-like domain structure for Ic 



Ic{Hcxt) — I I 



s vcl(I)h{Lx) sin(7r$a/$o) 



(21) 



Thus the dependence Ic{IIext) is given by the usual Fraunhofer curve with an effective critical current 
/co sin(7r$a/$o)/(7''^a/^o) the sign and value of which depends on dp, a etc. 
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Although in the theoretical papers [10-12] only a regular domain structure is considered, in real samples the domain 
structure is not strictly periodic. It may be almost regular (see for example, [21,22]), or very irregular ( [23,24]) with 
in-plane or out-of plane magnetizations. We study the effect of a possible irregularity of the domain structure on the 
basis of a simple model. We assume simply that the domain size a fluctuates around a mean value a and fluctuations 
of a are described by the Gaussian distribution. Then the dependence of /c(<5) on the dispersion of the fluctuations is 
given by the integral (in the absence of i?ea;t) 



where Ci = (/q°° daexp(— (a— a)^/(5^))~^ is a normalization constant. In Fig. 2 we plot the dependence /c(7r<l>a(a)/$o) 
for different parameter 7 = a/5. One can see that for large 7 this dependence coincides with a Fraunhofer pattern, 
but with decreasing 7 the amplitude of oscillations of Ic decreases and finally the function /c(^'o) does not change 
sign (no tt— states). 



In conclusion, using a simple model of a multidomain SFS Josephson junction, we have calculated the critical 
current Ic- It turns out that the current Ic changes signs when the in-plane magnetic flux = 4TrMa{2dF) in each 
domain equals n$o- The magnetic flux is caused by the magnetization in the ferromagnetic domains and therefore 
exists even in the absence of an external magnetic field. The oscillations of Ic observed experimentally by varying 
thickness 2dF or temperature T may be related not to the sign revcrsl of the local critical current density jc, but 
to a simple mechanism - the Fraunhofer-like oscillations of Ic caused by the internal magnetization M in domains. 
Almost nothing is known about the domain structure in real SFS junctions. For estimations we take AnMo « IkOe, 
a w Imkm, 2dp « 100^. For these values we obtain « 10~^Oe • crn^. This means that the critical current Ic 
changes sign for the thickness 2dp about 100^. This value of thickness is close to that used in experiments, although, 
strictly speaking, the values of the magnetization M and of the domain structure period a are not known. In order 
to make more convincing conclusions about what is the mechanism of the sign reversal effect (whether it is caused by 
the sign reversal of the critical current density jc or by the spatial phase modulation in a multidomain SFS structure) , 
further theoretical and, especially, experimental studies are needed. In particular, it would be interesting to study 
the influence of the domain structure on Shapiro steps in SFS junctions measured in a recent paper [20] 
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